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ABSTRACT 

Galactic interacting double white dwarfs (DWD) are guaranteed gravitational 
wave (GW) sources for the GW detector LISA, with more than 10 4 binaries 
expected to be detected over the mission's lifetime. Part of this population is 
expected to be eccentric, and here we investigate the potential for constraining 
the white dwarf (WD) properties through apsidal precession in these binaries. 
We analyze the tidal, rotational, and general relativistic contributions to apsidal 
precession by using detailed He WD models, where the evolution of the star's 
interior is followed throughout the cooling phase. In agreement with previous 
studies of zero-temperature WDs, we find that apsidal precession in eccentric 
DWDs can lead to a detectable shift in the emitted GW signal when binaries 
with cool (old) components are considered. This shift increases significantly for 
hot (young) WDs. We find that apsidal motion in hot (cool) DWDs is dominated 
by tides at orbital frequencies above > 10~ 4 Hz (10~ 3 Hz). The analysis of apsidal 
precession in these sources while ignoring the tidal component would lead to an 
extreme bias in the mass determination, and could lead us to misidentify WDs 
as neutron stars or black holes. We use the detailed WD models to show that for 
older, cold WDs, there is a unique relationship that ties the radius and apsidal 
precession constant to the WD masses, therefore allowing tides to be used as a 
tool to constrain the source masses. 

Subject headings: (stars:) binaries: general - stars: interiors - stars: white dwarfs 
- gravitational waves 



1. Introduction 



Binary systems containing compact objects constitute an important physical laboratory 
for investigating the properties of matter under extreme conditions. While direct inference 
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of these properties from "traditional" electromagnetic radiation is challenging or impossible, 
our knowledge about compact object properties can be greatly enhanced by measuring their 
interactions with their binary companions. The detection of gravitational waves (GWs) from 
these objects is expected to considerably advance our understanding of astrophysical systems 
harboring compact objects. Several ground-based GW observatories have already taken 
data a t design sensitivity (e.g., the Laser Interferometer Gravitational- Wave Observatory, 
LIGO; lAbbott et al.l 120071 ) . while the future generation of GW detectors inclu des the space- 
based Laser Interferometer Space Antenna (LISA; Danzmann et al. 1996, iHughed 12006 
and references therein). Our focus here lies on LISA, which was originally designed to 
be sensitive to GW frequencies between 10 -4 and 0.1 Hz, and to carry out high-precision 
astronomy for sources ranging from sub-solar-mass WDs to million-solar-mass black holes 
(BH). In its original design, LISA would have surveyed the whole galactic population of 
DWDs with period shorter than 6 h, with the potential of detecting f ro m 1,000 to more than 



10,000 sources ove r the lifetime of the mission ( iNelemans et al.ll2001aM l2004t iLiu et al.ll2010 



Ruiter et al.l 120101 ). providing the largest observational sample of these objects. At present 
it is not clear how the LISA sensitivity will be modified as a number of re-design options 
are being considered by the European Space Agency (ESA) and potential other partners. 
Regardless of the specifics it appears quite probable though that the sensitivity to Galactic 
DWDs at the middle to high-end frequency range relevant to this study will not be greatly 
affected. 

Most of the studies in the literature concerning DWDs have been focused on circular 
binaries, which are expe cted to domina t e the galac tic DWD populati on. However, recent the- 
oretical calculations by IWillems et al.l (120071 ) and iThompsonl (120101 ) predicted a population 
of eccentric DWDs, which could provide a unique opportunity for investigating degenerate 
matter with LISA. In fact, the properties of WDs can be studied t hrough the eff e cts o f 
tidal interactions on the detected GW signal from eccentric systems. IWillems et al.l (120071 ) 
studied eccentric binaries forme d through dynami cal interactions in globular clusters, (with 
eccentricities up to 0.8), while IThompsonl (120101 ) studied eccentric DWDs as products of 
hierarchical triple systems, where the tertiary induces Kozai osci llations in th e inner binary, 
driving it to high eccentricity (with eccentricities up to 1; see also lGouldll201ll for a discussion 
of the implications of such a population). 

In binaries with eccentric orbits, tidal forces cause a non-dissipative precession of the 
periastron of the orbit, also known as apsidal precession. This periodic motion leaves a mark 
in the GW signal that is potentially detectable. While in circular binaries GW radiation is 
emitted at multiples n of the orbital frequency z/, in the presence of periastron precession 
each of these harmonics is split into a triplet with frequencies nv ± ^{j-n and nu, where 
7 is the apsidal precession rate (IWillems et al.l 12008k hereafter WVK08). Various physical 
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processes can lead to apsidal precession; in the specific case of tides it occurs because the tidal 
deformation of the components perturbs the stars' external gravitational field, and therefore 
the Keplerian motion of the binary components. Tidally-induced apsidal precession depends 
on the binary orbital parameters and on the internal structure of the tidally deformed stars. 
Therefore, comparisons between observed and theoretically predicted rates can be used to 
place constraints on the WD properties. An equivalent approach has been very successful in 



the study of binary main -sequence stars (IClaret fc Gimeneall993l . I201CH ; I Claret fc Willems 



20021 ; ISchwarzschildlll958f ). 



Earlier studies carried out by WVK08 used zero-temperature WDs to investigate the 
contribution to apsidal precession rates due to general relativity (GR), tides, and rotation. 
WVK08 showed that, depending on the system parameters, tides could be the dominant 
contribution to apsidal precession in most of the LISA band. The aim of thi s paper is to 



expan d the initial study started by WVK08 and use detailed WD models from iDeloye et al. 



(120071 ). where the thermal properties of the WD are followed throughout the cooling phase, 
without zero-temperature assumptions. Even though we mainly focus our attention on 
detached DWDs, our results can be applied also to binaries containing a neutron star (NS) 
and a WD, which provide an even better probe to study WD physics. In DWDs apsidal 
precession is due to both components, while in NS-WD binaries the distortion of the NS 
due to tides and rotation contributes negligibly, and apsidal precession carries the unique 
signature of the WD. We devote a section of the paper to a discussion of NS-WDs. 

The plan of the paper is as follows. In §2] we summarize the equations governing the 
tidal, rotational and GR contribution to apsidal precession. In §3] we present the precession 
rates calculated using detailed WD models. We confirm the fundamental result previously 
obtained by WVK08 for zero-temperature models, and discuss the effect of WD temperature. 
In fjHwe study in more detail the tidally-induced apsidal precession rate. In §5]we discuss the 
implications of not properly accounting for tides when inferring WD masses from detected 
apsidal precession rates, while in §6] we investigate the source parameters that could be 
constrained from a detection. In SJ7] we discuss apsidal precession in NS-WD binaries. We 
conclude in SJH1 



2. Apsidal Precession In Binaries 

Here we summarize the equations governing the tidal, rotational, and GR contributions 
to apsidal precession in eccentric binaries. We consider a binary system consisting of two 
stars with masses Mi^, and radii -Ri,2, uniformly rotating with angular velocities Q\ y 2- We 
take the axes of rotation to be orthogonal to the orbital plane. Let 7 be the argument of the 



-4- 



periastron, P and a the orbital period and semi-major axis, respectively, and e the orbital 
eccentricity. Following WVK08, we consider orbital and rotational periods long com pared to 
the free harmonic periods of the component st ars (the so-called static-tides regime; ICowling 



59; 



19381 ; ISternd Il939l ; ISmeyers fa Willemsl 120011 ) . In this regime we can neglect resonances 



between the tidal forcing angular frequencies and the eigenfrequencies of WD free oscillation 
modes. Under these approximations, t he contributio n to the apsidal precession rate due to 
the quadrupole tides raised in star i is (jSternelll939l ): 



jTid,i 



30tt fRfY M- 



p 



Mi 



.1+32, 14 



1-e 



2^5 



(1) 



We assume that the WD components are sufficiently far apart that this quadrupole term 
dominates the precession (in §3.21 we show that higher order terms indeed contribute negligi- 
bly). As anticipated in JTJ 7™ depends on the orbital parameters, component masses, and 
internal structure of the tidally deformed star. The dependence on the internal structure is 
embedded in the quadrupolar apsidal precession constant fcj, which measures the star's cen- 
tral concentration. The constant ki ranges from for a point mass to 0.74 for an equilibrium 
sphere with uniform mass density. For any realistic stellar model ki assumes values between 
and 0.74. In standard notation the quadrupole apsidal precession constant is normally 
referred to as &2> where the subscript to k refers to the dominant longitudinal mode in the 
spherical harmonic expansion of the tide-generating potential (as explained in Sectj3j). To 
avoid confusion, we will drop the subscript 2, and use ki 2) to denote the apsidal 

precession constant for each binary component. 

A perturbation of the gravitational field and hence apsidal precession also arises because 
of the rotational distortion of the st ar caused by the centrifugal force. The rotationally- 
induced apsidal precession rate (e.g. ISternd Il939l ) depends on the star's internal structure 
and radius in a similar fashion as 7™: 



lR0ty 



2tt / Rj\ Mi + M 2 (tti/nf 

a ) Mi (1 - e 2 ) 2 h 



P 



(2) 



where Q = 1i\jP is the mean motion. 



Finally, the GR contribution to apsidal precession at the leading quadrupole order is 
given by: 

' 2tt ^ 5/3 3G 2 / 3 (M 1 + M 2 ) 2 / 3 



wher e G and c are the gravitational constant and the speed of light, respectively (jLevi-Civita 
19371 ). Note that 7^^ is independent of the internal structure of the components, and it is 
proportional to the total mass of the binary. 
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The total apsidal precession rate (7) is the sum of the three contributions above con- 
sidering both WD components. If a detection of the apsidal precession rate in a GW signal 
were to be used to infer the source's total mass, it is evident that not accounting for tides 
and rotation could potentially lead to an overestimate of the system mass, and a possible 
misclassification of the GW sources ($5]). 



3. Apsidal Precession Rates 



To investigate the various contributions to apsidal precession for eccentric DWDs in the 
LISA ba nd, WVK0 8 cons idered zero-temperature WDs following a mass-radius relation as 
given by lNauenbergl (119721 ). and setting the apsidal precession constant ki to 0.1 (appropriate 
for the WD masses considered in their study). Here we use the detailed WD structure 
models from lDeloye et al.l (120071 ) to compute ki at different evolutionary phases during the 
WD lifetime. 



3.1. Apsidal Precession Constant ki 

The apsidal precession constant is given by 

2k t = - 1, (4) 

Ki 

where is the radial component of the tidal displacement of a mass element of star i. 
In the limiting case of rotational and orbital periods long compared to the WD dynamical 
timescales (the so-called quasi-static ti des regime), £.i t r(r) is a so lution to the homogeneous 



second-order differential equation (e.g. ISmeyers fc Willemsll200ll ) 

d% T (r) | 2 / 1 dg(r) | 1 



dr 2 \g( r ) dr r 

X M-!fi±^ r (r) = 0, (B) 
dr r l 

where I is the longitudinal mode in spherical harmonics, and g is the gravity. The quadrupole 
apsidal precession constant is obtained by setting 1 — 2. The desired solution to Eq. (jSJ) 
should be finite at r = and at r — Ri satisfy the boundary condition 

<Kuir ^ ' ^(JJ,) = e i>T (2l + l)c l>0>0 . (6) 



dr J „ Ri 



Here e^T — {Ri/<^) 3 {Ms-i/Mi) is a dimensionless parameter indicating the strength of the 
tidal force with respect to gravity at the star's equator, while q^o are Fourier coefficients 
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of degree I, depending on the star's radius as (Ri/a) l ~ 2 . The dominant contribution of the 
terms associated with the second-degree spherical harmonics for systems with components 
sufficiently far apart justifies restricting our analysis to the tides generated by the terms 
I = 2. 



We consider He WD models from lDeloye et al.l (120071 ) with masses ranging from 0.1 M ( 



to 0.3 M Q in steps of 0.05 M©. For each mass we calculate the apsidal precession constant 
ki at different evolutionary stages. The results are summarized in Table [TJ while in Figure 
[TJwe show the time evolution of the apsidal precession constant for all the available masses 
(left panel). In what follows, we refer to the highest and lowest effective temperature (T e g) 
models as hot, and cool, respectively. For a given mass, a hot (cool) model represents a 
young (old) WD. We recall that ki measures the distribution of mass inside the star, and 
that it ranges from for a point mass to 0.74 for a uniform density sphere. As shown in 
Figure HJ ki increases with time, indicating that the WD becomes less centrally condensed 
while cooling. This trend is confirmed by the time evolution of the radial profile of the star's 
density, which becomes flatter as the WD evolves (right panel). 



3.2. An Example of Apsidal Precession Rates 

In Figure [2] we show an example of tidal, rotational, and GR apsidal precession rates 
calculated with our models, as a function of the orbital frequency v = 1/P for different 
orbital eccentricities and WD masses of 0.3 M©. The rotational contributions are calculated 
assuming pseudo-synchronism, i.e. that each WD's rotational angular velocity is equal to the 
orbital angular velocity at periastron (a note on this assumption follows below). The masses 
were chosen in order to compare our calculated apsidal precession rates against WVK08's 
results for zero-temperature WDs, with which our conclusions in general agree. We note 



Table 1. Apsidal precession constants of WD models 



0.1 Mq 


0.15 Mq 
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T cff (K) 


k 


T aB (K) 


k 


T aS (K) 


ft 


T cft (K) 


k 


5,448 


0.014 


9,743 
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5,340 
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0.127 
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0.134 


3,271 


0.122 
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0.113 
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0.153 
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0.142 


2,639 


0.134 


3,351 


0.127 


4,323 


0.122 






1,989 


0.146 


2,422 


0.139 


2,871 


0.133 


3,585 


0.127 
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that apsidal precession is dominated by GR for systems with small orbital frequency, while 
tides and rotation become the dominant contributions towards higher frequencies. As a 
reference, for the tidal contribution we include also the rates associated with the I = 3 and 



I = 4 longitudinal modes (jSterndll939l ). We note that the total phase shift to the signal from 
these terms is about a factor of 100 to 1,000 times smaller than the I = 2 tidal contribution. 
If synchronization at periastron is assumed, tides always dominate rotation. The same is 
true for a sub- synchronous system, given the dependency of Eq. <^ on the WD's rotational 
angular frequency. 

Observational constraints derived both from traditional spectrosco py and asteroseis- 
mology add to the growing evidence that single WDs are slow rotators (IHeber et al.lll997l ; 



Koester et al. 1999 



Berger et al.l 120051 ; iKawalerl 12004 ) . The only DWD system for which 



there has been an estimate of the ro tatio nal velocity of the components is S PSS J125733. 63+542850. 5. 
For this binary, iMarsh et al.l ( 120101 ) and iKulkarni fc van Kerkwijkl ( 120101 ) determined the ro- 
tation rate of both WDs to be super- synchronous, assuming the observed broadening of the 
spectral lines were solely due to rotation. However both studies cannot exclude that the 
broad ening seen could be due, instea d, to strong magnetic fields in more slowly rotating 
stars (IKulkarni fc van Kerkwijkll2010l and Marsh 2010, private communication). Further ob- 
servations are needed in order to settle this issue. Lacking more elaborate evidence that 
super-synchronous rotation is the norm, in this work we assume that synchronization at 
periastron is an upper limit on DWD rotation. 

It is clear from Figure |2] that apsidal precession is significant in most of the LISA 
band for both hot and cool DWDs, and high and low orbital eccentricities. At v ~ 1 mHz, 
e = 0.01, and for cool DWDs the precession of the periastron is expected to shift the phase 
of a GW signal by 3tt-Att over a 1 yr mission, and is therefore detectable (WVK08). The shift 
increases by more than a factor of 10 for hot WDs, and, at the same frequency of ~ 1 mHz, 
tides are the leading mechanism (accounting for ~ 85% of the total 7). For higher orbital 
eccentricities, tides dominate 7 at shorter orbital frequencies leading to much higher rates. 
We also note that, given the dependence of 7™,i on Rf binaries hosting hot (and bigger) 
WDs have 7™ systematically higher than cool (and contracted) DWDs by a factor of ~ 20: 
at v ~ lmHz and e = 0.01, j Tid)i for hot (cool) DWDs is ~ 7300(390) deg/yr. 

We note that among the available models, a DWD hosting 0.3 M & components yields 
the lowest jrid,i and 7R t,i, and the highest 7g_r,i- In fact, the larger radii of the lower mass 
models act to increase the tidal and rotational contributions relative to GR. The relative 
GR contribution, instead, being dependent upon the total system mass, decreases with the 
components' masses. 
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4. Tidal Contributions to Apsidal Precession 

GR-induced apsidal precession, if detected in binaries with long orbital periods (where 
this contribution is dominant), can be used to constrain the total system mass (see Eq. [3]). 
Here we show that 7™,i can be used in a similar fashion. 

Given the dependence of Eq. (CQ) on several stellar and binary properties (P, e, a, M±, M 2 , 
and Ri), an apsidal precession detection and rate measurement cannot be used to place tight 
constraints on the binary properties. In fact, although the orbital period and eccentricity 
can be precisely determined from the frequency spectrum of the GW signal, measuring radii 
and masses is much more challenging. In what follows, we show that some of the unknowns 
determining 7™,i can be eliminated if one considers the time evolution of the term kiRf 
that enters Eq. (pQ). As Figure [3] shows, during most of the WD's lifetime k{Rf assumes a 
plateau value which is unique for each mass. Consequently Eq. (JT]) effectively depends only 
on the orbital period, eccentricity, and component masses (through Kepler's law a can be 
re-written as a function of P, Mi, and M 2 ). According to our models, the plateau values of 
kiR\ (filled squares in Figure [3]) depend on the WD mass as follows: 

kiRf = -0.632 + 0.370 • Mr 1 - 709 , (7) 

where is the mass in solar units, and kiRf is in 10 _10 -Rq (see Figure H]). 

Even though for most of the WD lifetime there is a unique relation between kiRf and 
Mi that we can exploit to simplify Eq. ([I]), we note that the uniqueness of kiRf for a given 
mass is broken for young WDs. The probability density functions (PDFs) for the values of 
kiR 5 assuming an observational time uniformly distributed between and 10 Gyr after WD 
formation (consistent with our current understanding of star formation rate in the Milky 
Way) are shown in Figure [5J The dotted vertical lines represent the plateau of fcji2f for 
each mass, which is used for the one-to-one correspondence to the WD mass. In Table [2] 
we summarize the percent of overlap for the calculated PDFs of k{R\ , which represents the 
probability of confusing a less massive old (cool) WD with a more massive young (hot) WD. 
We conclude that, for a given fcji?f, we can distinguish between two masses within 0.05 M Q 
at the 90-95% confidence level, without any information on the age of the WDs. 

Given the unique dependence of kiR\ on the mass for cool WDs, the tidal contribution 
to apsidal precession can now be used to place constraints on the component masses also for 
short period binaries, extending this possibility to the entire LISA frequency domain. The 
substitution of Eq. (J7J) in Eq. ([[]) makes 7™,i depend on only P, e, and Mj, and since P 
and e can be precisely determined from the GW signal a measurement of 7 yields constraints 
on a particular combination of the masses in the system, which we call the "apsidal mass 
function." 
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Fig. 1. — Left: Apsidal precession constants as a function of time for different He WD 
models. Right: Radial density profile of a He WD with mass M; = 0.3 M at different 
evolutionary stages. The dashed line represents the hot WD model and evolution in time 
proceeds downward. The flattening of the profile as the star cools is common to all WD 
masses. 



Table 2. Probability of overlap in the values of kiR\ between an old, lower- mass WD with 
a young higher-mass WD, assuming observations uniformly distributed in the WD age. 
This probability is calculated from the time evolution of kiR\ appearing in Figure [3J 
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Fig. 2. — Apsidal precession rates for cool (left) and hot (right) DWDs with masses Mi = 
My, = 0.3 Mq, and orbital eccentricity e =0.01 (top), and e =0.5 (bottom). The tidal and 
rotational rates include the contribution of both stars. The dark- and light-grey dashed lines 
represent the tidal contribution from the longitudinal modes I = 3 and I = 4, respectively. 
Given the dependency of 7™,i and 7_R Q t,i on the contraction of the stars with time 
is responsible for the decrease of these two contributions with the evolutionary stage of 
the WDs. In the frequency range considered, the cool DWD at e =0.5 begins Roche-lobe 
overflow (RLO) at v ~ 2.7 mHz (at e = 0.01 the onset of RLO is at a higher frequency than 
displayed), while the hot DWD at e =0.01 (0.5) begins RLO at v ~ 1.8(0.65) mHz. 
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Fig. 3. — Evolution of kiRf as a function of the WD's age. Filled squares denote plateau 
values. The maximum age considered for each model is 10 Gyr. 
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Fig. 4. — Calculated fojiif as a function of the WD mass. Filled squares represent the plateau 
values of ktRf. The solid line is the fit given by Eq. (j7|). 
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Fig. 5. — PDF for kiRf for each model mass assuming an observational time uniformly 
distributed between and 10 Gyr after WD formation. The vertical dotted line represents 
the plateau of fcj-Rf for each mass. Note the non-uniform x-axis. 



-14- 



5. Misclassification of the Sources 



Ignoring tidal effects when trying to constrain the apsidal mass function from a measured 
rate of apsidal precession would lead to a dramatic overestimate in the masses, as shown 
in Figure |6j The left panel depicts lines of constant total 7 (i.e. GR, tides, and rotation) 
for cool DWDs with orbital frequency v = lmHz, and eccentricity e = 0.01. The tidal 
and rotational rates include the contribution of both stars. The right plot shows, for these 
same apsidal precession rates, the masses that would be inferred solely accounting for GR. 
Unless the distance to the source is known (see §6]), the amplitude of the signal gives no 
constraint on the masses, which would therefore be overestimated by orders of magnitude 
when only accounting for the GR parts of the apsidal mass function. In this case, an 
eccentric WD binary would be misclassified as a binary BH or NS. Such a misclassification 
would bias the ratios of different populations of compact object binaries that hold essential 
information about binary formation mechanisms and stellar evolution, and would hamper a 
correct astrophysical interpretation of the source's nature and properties. 



6. Adding 7 to GW Data Analysis: Constraining the Source Properties 



The G W signal emitted by DWDs at GW frequencies /(= 2y) ~ a few mHz is exp ected 



to show an intrinsic frequency evolution fjWebbink fc Hanlll998l ; IWillems et all 120071 ). We 
now discuss how a detection of apsidal precession in such binaries can be used to constrain 
the individual masses of the WD components. 

The measurable signal for GW detectors is the strain amplitude of the wave: 

y/9(n,e) 



h(n, e) 

AO 5 / 3 

M 



1.0 x 10 



-21 



?? 



P 
lhr 



-2/3 



d 



lkpc 



The pre-factor y/g(n, e)/n gives the sc ale of the strain ampl itude produced by the nth 



orbital harmonic at eccentricity e (e.g., iNelemans et al.ll2001bl ). P is the orbital period, d 
is the distance to the source, and M = (MiM 2 ) 3 / 5 / (Mi + M 2 ) 1//5 is the so-called chirp 
mass. If radiation reaction causes the orbit to evolve (and therefore the intrinsic frequency 
of radiation to change), assuming this drift is GR-driven, and that tidal and/or m agnetic 
spin- orbit coupling do not contribute, the chirp mass can be determined from (e.g. ISchutz 



19961 ) 



/ = 5.8 x 10" 



M 



5/3 



f^Hzs 



(9) 
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The distance to the source can then be calculated via Eq. (jSJ). It is clear that an apsidal 
precession detection combined with a determination of the chirp mass yields the masses of the 
individual components: M. and Eq. (jTJ) or (j3J) for tidally- or GR-induced apsidal precession, 
and M. together with Eq. ([TJ + dSj) in the regime where tides and GR contribute equally 
are sufficient to determine the individual masses of the components. Given the currently 
observed sample of rotation rates in isolated WDs and DWDs, we can assume that GR and 
tides are the only mechanisms contributing to apsidal precession. 

In the event that a signal has measurable /, the problem of source misclassification 
discussed in $5]is less acute. The chirp mass implied by the measured / would be inconsistent 
with the large system mass inferred from the GR contribution to apsidal precession alone; 
this inconsistency can be resolved, and the proper component masses determined, by properly 
including the tidal contribution to apsidal precession as discussed above. 



Neutron Stars - White Dwarfs Binaries 



More convenient laboratories for investigating WD physics via apsidal motion are NS- 
WD binaries. As opposed to DWDs, where the tidal and rotational contribution of both 
components have to be taken into account, apsidal precession in NS-WDs carries solely the 
signature of the WD, because the tidal and rotational distortion of the NS do not contribute 
to the apsidal precessi on. Eccentric NS-WDs are observed in the field typically as radio 



Kaspi et all 12000 



pulsars orbiting WDs (ILorimer 
fPSR B2303+46, 



Stokes et al 



Bailes et al. 



2005) , and two such binaries have b een discovered to date 
198,4 Ivan Kerkwijk fc Kulkarnil Il99fll and PSR J11416545, 



2003h . This field population of eccentric NS-WDs is thought 



to have formed from the evolution of binaries with similar mass components. The origi- 
nally more massive primary transferred mass to its companion to become a WD, while the 
secondary accreted enough matter to evolve into a NS, kicked at supernova into an eccen- 
tric orbit (ITutukov &: Yungelsonlll993l ). Even though no such systems have been observed 
in the LISA frequency band, theoretical simulations predict a wide range of orbital peri- 
ods and eccentricities , and a subset of this population might be readily detectable by LISA 
(IKalogera et al.ll2005l and references therein). 



Another advantage of using NS-WDs to investigate WD physics via apsidal precession 
lies in the higher number of expected systems if compared with eccentric DWDs in globular 
clusters. In fact, the expected formation rate of galactic eccentric NS-WDs is between 10 
to 10,000 times the expected formation rate of eccentric DWDs in globular clusters (see 
Kalogera et al.ll2005l and references ther ein for a comparison between empirical and theoret- 



ical rate estimates, 



Willems et all 120071 ) 
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Figure [7] shows, as an example, the tidal, rotational, and GR apsidal precession rates as 
a function of the orbital frequency for NS-WD binaries with different orbital eccentricities. 
To make a comparison with Figure [2] we consider a 0.3 M WD (Mi), while assuming a 
typical mass of 1.4 M & for the NS (M 2 ). The rotational contribution of the WD is calculated 
assuming synchronization at periastron. Again, GR dominates apsidal precession in binaries 
with the longest orbital period, while tides become progressively more important as the fre- 
quency increases. When compared to the DWD binaries shown in Figure [2j the dependency 
of Eqs. ([I]), (0), and ([3]) on Moleads lower tidally- and rotationally- induced apsidal motions, 
and a higher GR contribution.!]. Consequently, for the same orbital configuration and WD 
masses, the predominance of tides in a NS-WD occurs at higher orbital frequencies than 
in a DWD. Nevertheless, the tidal contribution needs to be taken into account to properly 
interpret the observations. 

A more detailed analysis of apsidal precession in eccentric NS-WD binaries will be the 
subject of future work. However, here we note that these systems potentially constitute an 
important astrophysical laboratory for testing our models. In NS-WDs, pulsar timing mea- 
surements might yield the masses of the individual components independently of precession 
effects, and therefore allow a comparison between the apsidal precession prediction of models 
and measurements. 



8. Summary, Discussion and Conclusions 

In this paper we use detailed WD models to continue the investigation started by 
WVK08 on apsidal precession due to tides, rotation and GR in eccentric DWDs. We find 
that apsidal motion can lead to a significant shift in the emitted GW signal, the effect being 
stronger for binaries with hot WDs, and weaker for cool WDs. 

In general, GR is the dominant mechanism at small orbital frequencies, while tides 
and rotation become increasingly important as the frequency increases. Given our current 
understanding of WD rotation rates, we can safely assume that apsidal precession is solely 
tidally-induced at frequencies where GR doesn't dominate the precession. Our detailed 
models also show that 7™,i in binaries hosting young WDs is systematically higher than in 
binaries hosting old WDs by a factor of ~ 20, but note that WDs spend only a small fraction 
of their lifetime at these temperatures. 



*E Mi,P, and e are kept fixed, using Kepler's law -y Tidti oc M 2 (Mi + M 2 )~ 5/3 , ju^i oc (Mi + M 2 )~ 2/s , 
and j GR ^ oc (Mi + M 2 ) 2 / 3 
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We investigate the astrophysical information that can be extracted from measured ap- 
sidal precession rates. First, we find that the component masses that would be inferred from 
the GR contribution alone could be overestimated by orders of magnitude. 

Secondly we show that, as the WD ages, fcji2f asymptote to a value that depends only 
on the WD mass. This behavior allows us to simplify the form of Eq. (pQ), and to use ^T%d,% 
to place constraints on the apsidal mass function for short-period binaries. As noted by 
WVK08, the leading GR contribution in long-period DWDs yields the total system mass, 
therefore our investigation extends the possibility of measuring the apsidal mass function 
to the entire LISA frequency domain. Furthermore, for systems where the distance can 
be independently measured via electromagnetic observations, or in which radiation reaction 
causes the orbit to evolve, apsidal precession and chirp mass measurements can be used to 
determine the masses of the individual components. 

Apsidal precession in eccentric NS-WD systems carries more straightforward information 
about the WD mass. These systems are formed in the field at rates 10-10,000 higher than that 
of eccentric DWDs in globular clusters. In these systems, the tidal and rotational distortion 
of the NS component contribute negligibly to 7, and apsidal motion carries the unique 
signature of the WD. We show that in NS-WD binaries apsidal precession can also lead to a 
considerable shift in the GW signal in most of the LISA band. The electromagnetic detection 
of these binaries as a WD orbiting a radio pulsar, combined with a possible determination 
of the components masses via timing measurements, provide a testbed for validating the 
reliability of our WD models. 

In this paper we investigate the importance of tidally-driven apsidal precession in low- 
frequency GW sources. However, tidal effects are also being discussed in the literature for 
high-frequency sources (/ > 10 Hz), such as NS-NS binaries (DNS). DNS are expected 
sources for the next generation of ground-base d GW detectors (e.g. Advanced LIGO; 



Harry & the LIGO Scientific Collaboration! |2010| ). In such binaries, the tidal distortion of 



the NS components and its imprint on the emitted GW signal during the inspiral phase 
could potentially be used to ex tract information about the NS equation of state (EOS), and 



therefore the NS matter (e.g. iFlanagan Hindered 120081 ; iHinderer et al.ll2010l ). Here we 



note that in the static-tides regime the leading mechanism to apsidal precession in DNS is 
GR, and that the tidal and rotational contributions to 7 can not be used to place constraints 
on the stars' central concentration. To give a sense of the magnitude of the various con- 
tributions to 7 in DNS, we consider a binary hosting two 1.4 M & NS emitting GWs at a 
frequency of 10 Hz (P = 200 ms). To maximize the tidal and rotational contributions to 
apsidal precession, we consider th e bigge st NS radius given by one of the commonly used 



EOS (~ 14 Km). IHinderer et all (120101 ) estimated the apsidal precession constant to be 
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&2 ~0.1. For this configuration of masses and radii, we consider the highest eccentricity not 
to have Roche- lobe overflow at periastron according to lSepinsky et all (120071 ). e ~0.85. Fur- 
thermore, to maximize the rotational contribution ( recall that A fR n t i i oc (Qj/Q) 2 ), w e consider 
a rotational period of 20 ms for both components (IReisenegger Goldreichl Il992l predicted 
the g-modes period in NSs to be as low as 10 ms). In this configuration, jcr is ~ 10 3 7™ 
and ~ 10 4 jRot- Obviously, the tidal contribution increases at the orbit shrinks. However, 
also decreasing the orbital period to 50 ms, while keeping the rotational period the same, 
7gr remains two orders of magnitude bigger than the tidal contribution. 

We conclude that the inclusion of tides is vital to the proper analysis of detected apsidal 
motion rates of DWDs and NS-WD binaries in the LISA band. DWDs, the most numerous 
guaranteed sources for LISA, hold essential signatures of binary and stellar evolution mech- 
anisms. The possibility of detecting apsidal precession in these binaries opens up a new and 
exciting window on WD physics that could not be investigated otherwise. 
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Fig. 6. — Lines of constant apsidal precession rate calculated for cool DWDs with orbital 
frequency v = lmHz, and eccentricity e = 0.01. Left: 7 calculated for different combi- 
nations of components' masses (Mi and M 2 ) accounting for GR, rotation and tides. Right: 
components' masses that would be inferred if only GR was taken into account. The various 
lines correspond to 7 in deg/yr of ~ 22570 (black-solid), ~ 5700 (dotted), ~ 2200 (dashed), 
~ 1100 (dot-dashed), ~ 680 (grey-solid). 
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Fig. 7.— Same as Figure El but for a NS-WD binary hosting a 1.4 M NS and a 0.3 M 
WD. In the frequency range considered, the cool WD at e =0.5 begins RLO at v ~ 2.5mHz 
(at e = 0.01 the onset of RLO for a cool WD is at a higher frequency than displayed), while 
the hot WD at e =0.01 (0.5) begins RLO at v ^ 1.7(0.6) mHz. 



